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ABSTRACT 

We analyse the dependence of clustering properties of galaxies as a function of their large- 
scale environment. In order to characterize the environment on large scales, we use the cat- 
alogue of future virialized superstructures (FVS) by Luparello et al. and separate samples of 
luminous galaxies according to whether or not they belong to FVS. In order to avoid biases 
in the selection of galaxies, we have constructed different subsamples so that the distributions 
of luminosities and masses are comparable outside and within FVS. As expected, at large 
scales, there is a strong difference between the clustering of galaxies inside and outside FVS. 
However, this behaviour changes at scales r < 1 h Mpc, where the correlations have sim- 
ilar amplitudes. The amplitude of the two-halo term of the correlation function for objects 
inside FVS does not depend on their mass, but rather on that of the FVS. This is confirmed 
by comparing this amplitude with that expected from extended Press-Schechter fits. In order 
to compare these observational results with current models for structure formation, we have 
performed a similar analysis using a semi-analytic implementation in a ACDM cosmological 
model. We find that the cross-correlation functions from the mock catalogue depend on the 
large-scale structures in a similar way to the observations. 

From our analysis, we conclude that the clustering of galaxies within the typical viri- 
alized regions of groups, mainly depends on the halo mass, irrespective of the large-scale 
environment. 
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1 INTRODUCTION 

The large-scale structure of the Universe, and in particular the 
largest virialized systems, are a key probe of the evolution of 
the density field from the primordial fluctuations, and conse- 
quently of the nature of the energy content of the Universe and 
the long range action of gravity. Thanks to the widely available 
galaxy redshift surveys performed in the last years, as Las Cam- 
panas Redshift Survey l lShectman et al.ll 19961). the 2-de gree Field 



Galaxy Redshift Survey (2dFGRS 



Collcss et al 



Sloan Digital Sky Survey (SDSS. lYork et alj200d : 



2001) and the 



Stou ghton et al.l 



l2002t) . we know that the large-scale structure of the Universe ap- 
pears as a network made up of walls, filaments, knots and voids 
djoeve er et al. 1978; Gregory & Thom psor] 1 19781 : IZeldovich et al.l 
1 19821 ; Ide Lapparent et al.ll 19861) . The nodes are the intersections of 
walls and filaments, so they are the highest density regions. It is 
these regions, or some part of them, which we call superstructures. 
In a ACDM cosmological model, where the dynamics of the Uni- 
verse is expected to be dominated by the accelerated expansion, 
superstructures can be defined as the currently overd ense regions 
that will be bound and vir i alized systems in the fut ure tiBusha et al. 
120051 : iDtinner et al]|2006l ; lArava-Melo et alj|2009b iLuparello etaU 



l201ll) . These superstructures provide important information about 
the matter distribution on cosmologica l scales, allowing precise 
analysis of the cosmological model ( Kolokotronis et al] 20021 : 
lEinasto et al] l200l lArava-Melo et al] 120091 : lEinasto et alj [201TI ; 
Sheth & Diaferia feoill) . The clustering properties of galaxies on 
scales smaller than the size of superstructures are key to observa- 
tionally c onstrain the accretion processes that give rise to luminous 
galax ies JSpringel et~al]|2005b [Bildfell et aT]|201 lb iFontanot et all 
2009). In this context, the local environment is claimed to be 
the principal factor in defining the global properties of galaxies 
and their distribution dGriitzbauch et al. 201 lb Kimm et al] 1 20091 : 
iBlanton et al]|2005b|Park et al]2007b|Park & Choill2009l) . 

In fact, several models of galaxy formation assume that galaxy 
properties are determined by the haloes where they for med, and not 
by th e large scale environment that surrounds them jBerlind et al] 
120031 : lYangetal]|2003l : iBaugr] 12006b iKauffmann et al]|l997l) . In 
this context, the population of galaxies in a halo of a given mass 
should be independent of the location of the halo. However, in the 
last years there were different studies, both observational and using 
simulations, which show that the galaxy properties, such as the star 
formation rate or galaxy colours, depend on large-scale structure 
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I Park & Choil2009l;lBinggelill982l ; lDonoso et alj2006l : ICrain et al.l 



l2Q09l ; IWhiteetalj2010l) 



The formation and evolution of systems that are embedded in 
superstructures could be conditioned by these l arge overdensities 
drloffman et all 120071 ; lArava-Melo et all 120091 ; IBond et aljfcoid 
and references therein). Galaxies within them also show different 
properties and spati al distributions (lEinasto et alj|2003l;IWolf et al l 
20051; lHaines et alj l20od: lEinasto et al] l200l IPorter et alj 120081; 
Tempel et al. 2009 ; Fle gnor & Johnston-Holl itt 2010; Tempel et al.l 



201 II ; lEinasto etal.11201 II) . 



Recent galaxy redshift surveys J York et al.l200 0; Colless ^t al.1 
l200lh sample a sufficiently large volume to allow the study 
of galaxy formation and evolution from a statistical perspec- 
tive and unveil its details with increasing confi dence. The aim 
of th is work is to use SDSS-DR7 catalogue (Abazaii an et al.1 
2009) to study the influence of superstructures on galaxies. We 
will study the dependence of the clustering of faint galaxies 
(tracers, -21.0 < M r < -20.5) around brighter galaxies (centres, 
-23.0 < M r < -21.0) according to whether the centre galaxies are 
part of superstructures. In this work, we will analyse the clustering 
properties of galaxies located in large structures which are still un- 
dergoing a virialization process. Then we can study the differences 
in the small scale clustering of these galaxies compared to galaxies 
in the field. For this st udy, we use the superstructures identified by 
iLuparello et al. I feOllh . dubbed F uture Virialized St ructures (FVS), 
identified in the SDSS-DR7 an d lZapata etai] f2009) groups which 
will be used to characterize the virial mass of the groups which host 
our galaxies. 

This paper is organized as follows. In Section [2] we describe 
data collected for groups of galaxies and future virialized struc- 
tures. We then analyse the clustering properties using the cross- 
correlation function, as described in Section[3] and show results in 
Section|4] We summarize and discuss the results in Section|5] 



2 DATA AND SAMPLES 

2.1 SDSS-DR7 Galaxy Catalogue 

The Sloan Digital Sky Survey dYork et alj200oh is one of the largest 
and most ambitious surveys carried out so far. It has deep multi- 
colour images covering more than one quarter of the sky, and spec- 
tra for about 930000 galaxies and 120000 quasars. It is one of the 
largest data sets produced and contains images, image catalogues, 
spectr a and redshifts. In this work we use the Seventh Data Release 
(DR7. lAbazaiian et alj|2009h . which is publicly available^. 

The spectroscopic galaxy catalogue comprises a footprint area 
of 9380 sq.deg. The limiting apparent mag nitude for the sp ectro- 
scopic catalogue in the r-band is 17.77 dStrauss et al.ll2002l) . We 
use a more conservative limit of 17.5 to ensure completeness. Also, 
in order to avoid saturation effects in the photometric pipeline, we 
consider galaxies fainter than r = 14.5. The image deblending soft- 
ware often fragments images of bright galaxies with substructure, 
so our cuts prevent possible artifacts in the final catalogue^. 



2.2 SDSS-DR7 Superstructures 

The maps of the Univ erse depict a complex network of fil- 
aments and voids (e.g. lEinasto et aT] 1 19961; IColless et al] l200ll ; 



1 http://www.sdss.org/dr7 

2 http://www.sdss.org/dr7/products/general/target_quality.html 



|jaanisteetai]|2004 lEinastd I2OO6I ; I Abazaiian et"ai] [2009), where 
clusters of galaxies are preferentially located at the nodes of fila - 
mentary structures dGonzalez & Padillal20ldlMurphv et al.l201 ll) . 
This picture was possible only after extended galaxy surveys were 
completed, and it is also supported by numeric al experiments con- 
sistent with the standard cosmological mod el I Frisc het al.1 1 1 995l ; 
IBond et alJl996l ; ISheth et al]2003£lShandarin et alj2004h . 

Large scale structures, which we will refer to as superstruc- 
tures, have been studied with a variety of methods. The first statisti- 
cal studies of superstructures were performed by li nking Abell clus- 
ter positions dZucca et aljl9 93; Einas toet aljl997h . Later, the real- 
ization of wide-area sur veys of galaxies with r edshift as Las Cam- 
panas Redshift Survey dShectman et al.ll 19961). the 2-de gree Field 
Galaxy Redshift Survey (2dFG RS, IColless et al.l l200lh and the 



Sloan Digital Sky Survey (SDSS. lYork et alj2OO0l : 



Stoughton et al. 



2002), allowed the identification o f superstructures dir ectly from 
the large-scale galaxy distribution. lEinasto et alj ( 120071) identified 
superclusters in the 2dFG RS using a density field method, and 
ICosta-Duarte et all d201 ll) studied the morphology of superclusters 
of galaxies in the S PSS. The largest catalog ue of superclusters has 
been constructed bv lLiivamagi et al] d2012 ). who implemented the 
density field method on the SDSS-DR7 main and luminous red 
galaxy samples. In all cases, superclusters are operationally defined 
as objects within a region of positive galaxy density contrast, and 
thus are subject to a certain degree of arbitrariness in their detec- 
tion. 

Within the ACDM Concordance Cosmological Model, an ac- 
celerated expansion dominates the present and future dynamics of 
the universe and thus determines the nature of gravitationally bound 
structures. Therefore, an alternative definition of these large-scale 
structures is that of overdense regions in the present-day universe 
that will become bound and virialized structures in the future. Thus, 
under the assumption that the luminosity is a somewhat unbiased 
tracer of mass on large scales, the integrated luminosity density 
of galaxies is c ommo nly used as an indicator of mass density. 
Luparello et al.1 d201 ll) applies the luminosity density field method 
to identify FVS in large spectroscopic galaxy surveys . The id entifi- 
cation procedure is based on the theoretical criteria of Dtin ner et al] 
(2006). Using numerical simulations, they establish the minimum 
mass o verdensity necc e sary for a structure to remain bound in the 
future. lLuparelloetai1d201ll) use this mass overdensity criteria to 
calibrate the luminosity-density threshold, and identify superstruc- 
tures as systems that are likely to evolve into virialized structures 
in the distant future, in the assumed cosmology. The catalogue 
of FVS was extracted from a volume-limited sample of galaxies 
from the SDSS-DR7, in the redshift range 0.04 < z < 0.12, with 
a limiting absolute magnitude of M r < -20.47. The luminosity- 
density field is constructed on 1 fr 1 Mpc cubic cells grid, ap- 
plying an Epan e chnik ov kernel of r = 8 Ir 1 Mpc (equation 3 of 
ILuparello et al] d201lh ). The luminosity overdensity threshold is 
fixed at D T = pi, m ,/pi„„, = 5.5, and the structures are constructed 
by linking overdense cells with a simple version of a "Friends of 
Friends" algorithm. In order to avoid contamination from smaller 
systems, they also assign a lower limit for the total luminosity of 
a structure at L„,„ c , > 10 12 L o . The main catalogue of superstruc- 
tures has completness over 90 per cent and contamination below 5 
per cent, according to calibrations made using mock catalogues. 
The volume covered by the catalogue is 3.17 x 10 7 (fr 1 Mpc) 3 , 
within which 150 superstructures were identified, composed by a 
total of 11394 galaxies. FVS luminosites vary between 10 12 L Q 
and - 10 14 L Q , and their volumes range between 10 2 (IT 1 Mpc) 3 
and 10 5 (fr 1 Mpc) 3 . Because of the luminosity density field depen- 
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dence, Luparel lo et al. l|) analysed 3 samples with different lu- 
minosity thresholds, dubbed SI, S2 and S3, described in table 1 
of their paper. We will consider sample S2 in our analysis, which 
contains 89513 galaxies with M r < -20.47 in the redshift range 
0.04<z<0.12. 

2.3 SDSS-DR7 Galaxy Groups 

The aim of this paper is to analyse the clustering of galaxies con- 
sidering different samples. Some of these samples are selected tak- 
ing into account the mass of the systems considered. To estimate 
the mass distributions of the samples we use the estimated mass 
of galaxy groups. The galaxy groups correspo nd to those identifie d 
in the SDSS galaxy catalogue pres ented by [ Zapata et alj d2009b . 
extended to cover the SDSS-DR7 . IZapata et alj d2009b identifie d 
groups using the same method as Merchan & Zandivarezl d2005t) . 
with a "Friends of Friends" algorithm, with a variable projected 
linking length and a fixed radial linking length. These were set by 
Merchan & Zandivarez to obtain a sample as complete as possible 
and with low contamination (95 per cent and 8 per cent, respec- 
tively). The catalogue contains 83784 groups with at least 4 mem- 
bers, and is limited to redshift z < 0.12. 



2.4 Synthetic catalogue and simulation 

In order to test the procedures implemented in this work, we 
constructed a mock galaxy catalog ue from a semi-ana lytic model 
of galaxy formation (GALF ORM, iBower et al.1 2008) applied to 
the Millennium simulation dSpringel et alj 120051) , which adopts 
a ACDM concordance cosmological model. The semi-analytic 
model of galaxy formation collects information from the merger 
trees extracted from the simulation, and generates a population of 
galaxies within the simulation box. The mock catalogue is set up 
to mimic the geometry of the SDSS-DR7 footprint and reproduces 
the dilution in the number of galaxies with redshift. 



3 METHOD 

The clustering of galaxies is analysed using the two-point corre- 
lation function. Early studies of the spa tial distribution of galaxies 
using this technique were presented by IPeeblesI dl980l). The two- 
point correlation function, £(r), is defined as the excess of the prob- 
ability, SP , of finding a galaxy, defined as tracer, at a given distance 
from another galaxy, dubbed the centre. Thus, 



6P = nM + %(r))dV, 



(1) 



where dV is a volume element and n g is the mean number den- 
sity of tracer galaxies . There are several estimators of this function 
dKerscher et alfe OQO ) that compute this probability by counting ob- 
ject pairs. In this work, we use one of the most commonly used es- 
timators to determine th e cross-correlation function, as defined by 
|Pavis& Peebles! dl983b : 



Pi£> 2 N Rl N R2 
RiR 2 N D .N Dl 



1. 



(2) 



where DjDj is the average number of object pairs and R t Rj is the 
average number of pairs in a random sample. N D . and N R . represent 
the number of objects in the data catalogue and in the random sam- 
ple, respectively. The random sample is generated with the same 
geometry as the real catalogue using the same angular and radial 
selection functions. In particular, the advantage of this estimator 



is that if the sample D l has few objects, R l can be selected to be 
larger than Di, and thus minimize the noise in the <f calculation. 
Al though there are more accurate estimators, like the one defined 
bv lLandv & Szalavl d 1993b . differences between both esti mators are 
negligible due to the large volume size of the SDSS (e.g. lPazetall 
1201 lb . 

In galaxy redshift catalogues, where the distance is estimated 
using the spectroscopic redshift which includes a peculiar veloc- 
ity component, the three dimensional distribution of galaxies is af- 
fected by distortions in the line of sight. To take this into account 
we estimate the correlation function £(cr,n), as a function of the 
projected (cr) and line of sight (n) distances. As these distortions in 
redshift space occur on ly in the radial directi on, we implement the 
inversion presented by ISaunders et alj d 1992b to obtain the spatial 
correlation function £;(r) from £(cr, re). We integrate along the line 
of sight to obtain the projected correlation function H(cr): 



H(cr): 



£(<r, n)dn = 2 



Jo 



cr 2 + y 2 )dy. 



(3) 



In this work we use cr = 0.3 and x = 0.3 Mpc as lower integration 
limits for this equation, and cr = 30 and n = 30 Mpc as upper limits. 
We use diferent binning schemes for £(cr, jx) according to the size 
of each sample, being the number of bins in the range 13 to 19. 

Then, we can directly estimate the real space correlation func- 
tion (£(>)) by the inversion of E(cr) assuming a step function 
H(cr,) = Hj in bins cen tered in cr, and interp olating between r = a, 
values (equation 26 of ISaunders et al.ll 1992b : 



0>i - crj 



In 



'7+1 



7+1 



(4) 



This equation is a simple and direct way to invert the projected 
correlation function. 

In addition, we also compute the correlation function in red- 



shift space, £;(s), where s 



• + y — IxycosO is the redshift-space 



distance; and x and y are the line of sight distances to each object 
of the pair, with angular separation given by 6 = cr/y. £(s) con- 
tains information about the peculiar velocities of the galaxies in the 
sample. 



4 RESULTS 

We sele ct galaxies which are m embers of the superstructures de- 
fined by Lupare llo et alj d201 lh to construct the samples of galax- 
ies in FVS, while we select galaxies which are members of the 
SDSS-DR7 galaxy catalogue but are not members of FVS for the 
samples of galaxies outside superstructures. To analyse the cluster- 
ing properties of both samples, we measure the cross-correlation 
function, using faint galaxies with luminosity -21.0 < M r < -20.5 
as tracers. In both cases we select galaxies with luminosity 
-23.0 < M r < -21.0 as centres. The centre and tracer galaxies are 
above the luminosity limit for the volume limited sample of the 
SDSS-DR7 galaxy catalogue. 

In Section |4~T1 we analyse the correlation functions around 
centre galaxies inside and outside FVS. The difference in the clus- 
tering of faint galaxies around bright galaxies is not only produced 
by the selection of galaxies according to their FVS membership, 
but also is related to other processes that occur on smaller scales. In 
particular, the luminosity of galaxi es is known to correlate with the 
clustering amplitude of ga l axies |Alimi et al] 



l2005t ISwansonetal] |2008: Wa ng et alJl2oTT 



1988; Zehavi et al. 



Zehavietalj|2011l : 
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Mont 


-23.0 < M r < -21.0 
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5563 


33306 


-23.0 < M r < -21.0 




8G 




2013 


4603 


-23.0 < M r < -21.0, same luminosity dist. 


gL 






5102 


30526 


-23.0 < M r < -21.0, same mass dist. 








1765 


3900 


-23.0 < M r < -21.0, same mass and luminosity dists. 




g G ML 




1765 


3316 


-23.0 < M r < -21.0, same mass and luminosity dists., in mock 




go ML — mock 




5550 


8711 


> 8 members and same mass distribution 

> 8 members and same mass dist., in mock 






GSM 
GSM - mock 


509 
708 


1829 
2145 



Table 1. Description of samples of galaxies or groups used as centres in the computation of f(r). Each name indicates objects (galaxies or groups) inside 
and outside FVS by the "in" and "out" sufixes, respectively. When indicated, subsamples are selected so that the distributions of galaxy luminosities or group 
masses inside and outside FVS are comparable. The last two columns indidate the number of objects in each centre sample. 



iRoss et al J201 ll> . There is also evi dence of the depend ence of clus- 
tering with dark matter host mass feehavi et al]2012h . Then, in or- 
der to disentangle and quantify the contribution of luminosity and 
mass to the total correlation amplitude, we will analyse the effect 
of equating luminosity and mass distributions on the correlation 
signal in Sections |4.2| and |431 In order to perform this analysis, 
we will use different samples of centre galaxies selected according 
to different properties. In Table Q] we give a brief summary of the 
samples (galaxies or groups) which are used as centres in the com- 
putation of correlation functions. Each sample name is used for two 
samples: one for galaxies located in a FVS, which we denote with 
the "in" subindex, and the other for galaxies which are not part of 
any FVS, identified with the "out" subindex. The number of objects 
in each sample is also given in TableQ] 

In Section |4~Tl we will analyse samples g,„ and g out , which con- 
tain bright galaxies (-23.0 < M r < -21.0) of the SDSS catalogue. 
In Section l4~2l we will analyse samples gL which are defined from 
the g samples, but with the additional restriction that their lumi- 
nosity distributions are comparable. In Section l4~3l we introduce 
the SDSS-DR7 galaxy groups in the analysis and we separate the 
groups located inside and outside FVS. Then, we define the sam- 
ples ga selecting bright galaxies contained in these groups. 

To analyse effects on the galaxy selection, in Section |4~4l we 
define new samples from the group catalogue. We redefine the 
group samples of Section 14.31 in order to have the same group 
virial mass distributions. The samples denoted by g G M are obtained 
by selecting the bright galaxies contained in these groups. To fur- 
ther constrain our samples, in Section POl we redefine the samples 
go ML by adding the condition that the bright galaxies located in 
groups with the same mass distributions also have the same lumi- 
nosity distributions, i.e, the gGML samples contain bright galax- 
ies with the same luminosity distributions located in groups with 
the same mass distributions. In Section |4~6l we select galaxy groups 
with at least 8 members, and then divide them according to whether 
they are located in FVS or not, and we make their mass distributions 
similar by selectively limiting the sample. In this case, to define the 
G8M samples we select geometrical centres of groups instead of 
the bright galaxies, so that we can use virial mass estimated for 
groups to uncover mass effects. In Section |4~7l we define the same 
samples as gcML and G8M but taken from a mock catalogue. We 
name these samples as g G ML-mock and GSM -mock respectively, 
and use them to test the reproducibility of our results by current 
models for structure formation. 
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Figure 1. Cross-correlation functions of SDSS-DR7 galaxies for the g 
samples. The dashed lines correspond to galaxies in sample g,„ and the 
solid lines correspond to galaxies in sample g out . 



4.1 Clustering of faint galaxies around bright galaxies 

Fig. [T] shows the cross-correlation functions for samples g,„ and 
g mll , i.e., all galaxies in the luminosity range -23.0 < M r < -21.0. 
The dashed curves show the results obtained for the sample of 
galaxies in FVS, g,„, while the solid lines show the correlation func- 
tion obtained for the sample of galaxies not contained in FVS. The 
shadowed regions indicate cosmic variance estimates computed us- 
ing Jackknife statistics over the FVS catalogue. The number of 
Jackknife realizations is equal to the number of FVS in each case. 
As can be seen from this plot, the probability excess of finding a 
centre-tracer pair of galaxies is higher for the sample of galaxies 
contained in FVS than that corresponding to the sample of galax- 
ies outside FVS, i.e. the clustering of galaxies is greater when they 
are contained in superstructures. A difference in clustering on large 
scales (two-halo term) is expected given the selection criteria ap- 
plied to the samples. Galaxies in the same range of brightness are 
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Figure 2. (a) Luminosity distributions. The left panel corresponds to the g samples and the right panel shows the luminosity distributions of the new samples 
(gL), defined to have the same luminosity distributions. In both cases, the dashed curves correspond to the centres galaxies belonging to FVS while the solid 
curves correspond to the centre galaxies outside FVS. The bottom panels show the empirical cumulative luminosity distributions for each case, (b) Cross- 
correlation functions of SDSS-DR7 galaxies for gL samples. The dashed lines correspond to galaxies in FVS and the solid lines correspond to galaxies outside 
FVS. Both samples have the same luminosity distributions. 



more clustered if they are in FVS, on small and large scales (one- 
halo term and two-halo term). 

From Fig. Q] it is evident that the clustering of galaxies is 
greater for galaxies located within the FVS, which suggests the pos- 
sibility that the formation histories of galaxies and their collapse, 
are strongly influenced by the large-scale environment. However, 
this signal should be taken as an upper limit for this effect. Indeed, 
other correlations caused by differences in the galaxy selection of 
the two samples could diminish the clustering difference. In order 
to quantify how this selection affects the results, we will consider a 
number of restrictions on galaxy properties, using other samples of 
TableQ 



4.2 Clustering of faint galaxies around bright galaxies: 
uncovering galaxy luminosity effects 

The clusteri ng of galaxies is known to depend on the luminosity 
of ga l axies dAlimi et al J 1 19881: IZehavi et alj|2005h ISwanson et al.l 
l2008l : IWang et alj 1201 ll ; IZehavi et alj l201ll ; iRoss et alj l201ll) . 
Therefore we will analyse how this dependence affects the previous 
results. The median luminosity of g oul galaxies is fainter than that of 
gi„ galaxies. This means that centres outside FVS are typically less 
luminous, and therefore their clustering amplitude would be lower 
than that of galaxies which are located in FVS due to this. To rule 
this out, we define samples gL in and gL m „ from the original g sam- 
ples, trimming objects so that the final luminosity distributions are 
similar. This proccess is made by comparing the original luminos- 
ity distributions and then randomly removing objects from either of 
these two samples. This removal is proportional to the differences 
between the distributions at a given luminosity. We progressively 
discard objects until the two distributions are reasonably similar. 



The left panel of Fig. |2(a)| shows the luminosity distributions 
of the original g samples, where dashed curves correspond to sam- 
ple gj„ and solid curves correspond to the sample g mt . In the right 
panel we show the luminosity distributions of the new samples gL, 
trimmed to have the same luminosity distributions, where dashed 
curves correspond to gL in and solid curves correspond to gL m „. The 
bottom panels show the empirical cumulative luminosity distribu- 
tions for each case. 

We estimate the cross-correlation functions for both gL sam- 
ples, shown in Fig |2(b)| Here we use the same line style coding of 
FigureQ] As can be seen there is a statistically significant difference 
between the clustering amplitudes of both samples, being greater 
for galaxies contained within FVS. This signal is observed at all 
scales including small scales (one-halo term), suggesting an effect 
of the surrounding FVS on the galaxy environment. This signal is, 
by construction of the samples gL, independent of the luminosity 
of centre galaxies but it still may depend on the mass of haloes they 
inhabit. 



4.3 Clustering of faint galaxies around bright galaxies: 
Galaxies in groups 

In this section we study the origin of the previous signal of dif- 
ferent clustering amplitude for galaxies of equal luminosity inside 
and outside FVS, focusing on the virial mass of the host groups 
of galaxies, which can be interpreted as probing larger, but still in- 
termediate scale structures. We select galaxy groups for which we 
have estimates of their mass. As a first step, we compare the cross- 
correlation function of galaxies in groups, distinguishing between 
groups that are located inside and outside FVS. To this end we use 
the SDSS-DR7 galaxy group catalogue described in Section 1231 
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r [Mpc] 

Figure 3. Cross-correlation functions of galaxies in groups (go samples). 
The dashed lines correspond to galaxies in ga in sample and the solid lines 
correspond to galaxies in ga„„, sample. 



We define two new samples of centre galaxies contained in these 
groups, one with groups belonging to the FVS (go Ul ) and the other 
one with groups that do not belong to any FVS (gc a „,)- Galaxies 
are again restricted to luminosities within -23.0 < M, < -21.0, as 
in the previous analyses. For a more direct comparison, we use the 
same tracer sample as in previous sections. 

The cross-correlation functions of the centre-tracer pairs for 
both samples are shown in Fig. [3] The dashed curves correspond to 
sample g Gjn , while the solid curves correspond to galaxies in sample 
gG mr 

4.4 Clustering of faint galaxies around bright galaxies: 
uncovering dark matter mass effects 

In this section, we analyse the possible effects coming from differ- 
ences in the mass of the groups where the centre galaxies reside. 
Therefore, we repeat the analysis performed for the luminosity de- 
pendence (Section \4.2\ , redefining instead samples with the same 
mass distributions. To determine the mass distributi ons we use 
the SD SS-DR7 galaxy groups catalogue described in Zap ata et al.l 
(2009) as in the previous section. We define two samples of galaxy 
groups, one with groups belonging to FVS and other one with 
groups outside FVS. The median of the mass distribution of the 
galaxy groups that do not belong to FVS corresponds to a lower 
value than the median of the mass distribution of groups inside 
FVS. We adjust these samples and redefine two new samples of 
groups, trimmed so that their mass distributions are similar. 

The left panel of Fig. [4(a)] shows the mass distributions for the 
original samples, where dashed lines correspond to groups located 
in FVS and solid lines correspond to groups outside FVS. With the 
same line style coding, the right panel shows the mass distributions 
of the new samples, redefined to have the same mass distributions. 
The bottom panels show the empirical cumulative mass distribu- 
tions for each case. 



From the latter two samples of groups, we define two samples 
of centre galaxies contained in these groups: one with galaxies in 
groups belonging to FVS (goM,„) and the other one with galaxies in 
groups outside FVS (gcM oul ). Fig. |4(b)| shows the cross-correlation 
functions of these samples. The dashed lines correspond to the sam- 
ple gaM in while the solid lines correspond to the sample gc,M om . 
As can be seen from this figure, the correlation function amplitudes 
differ only slightly in the inner regions (r < 1 Mpc). This is a similar 
behaviour to that observed in the case of centre galaxies in groups 
(-23.0 < M r < -21.0), selected without restrictions in group mass. 



4.5 Clustering of faint galaxies around bright galaxies: 
combined effects from luminosity and mass 

In this section we analyse how the luminosity and mass dependence 
of clustering affect the previous results. To do this, we repeat the 
previous process restricting samples to have comparable galaxy lu- 
minosity and host group mass distributions simultaneously. 

The left panel of Fig |5(a)| shows the luminosity distributions 
of the centre galaxies of the samples g G M described in previous 
section, which are defined to have similar mass distributions. The 
dashed curves correspond to galaxies in the sample gc;Mj„ while 
the solid curves correspond to galaxies in the sample gaM mll . The 
middle panel of Fig. |5(a)| shows the luminosity distributions of the 
new samples (go ML,,, and g G ML m „) defined by equating the lumi- 
nosity distributions of the previous samples (same line coding). The 
bottom panels show the empirical cumulative distributions for each 
case. 

It is important to ensure that the mass distributions of the hosts 
of centre galaxies remain similar once the luminosity distributions 
are trimmed to be comparable. The right panel of Fig |5(a)| shows 
the mass distributions of the groups of the samples g G ML. The 
line coding is the same as in the previous plots. From this plot 
we can corroborate that the mass distributions of both samples re- 
main comparable. The bottom panel shows the empirical cumula- 
tive mass distributions for each case. Thus, we obtain two sam- 
ples of centre galaxies with similar luminosity distributions, pop- 
ulating groups with comparable mass distributions. One of these 
samples contains galaxies in groups which are members of FVS 
(gcMLj„), and the other comprises galaxies in groups away from 
FVS (gcML mll ). The cross-correlation functions of the centre- 
tracer pairs for both samples are shown in Fig |5(b)| The dashed 
curves correspond to galaxies in sample go ML,,,, while the solid 
curves correspond to the galaxies in sample g G ML mll . 

From this figure we conclude that at small scales, the ampli- 
tudes of the clustering of both samples are slightly different, but 
within a l-cr uncertainty. 



4.6 Clustering of faint galaxies around centres of groups: 
GSM samples. 

In this section we analyse the group-galaxy cross-correlation func- 
tion, considering faint galaxies as tracers, as in the previous cases, 
but in this case selecting geometrical centres of groups instead of 
the brightest group galaxies, i.e, we study the dependence of the 
clustering of faint galaxies around centres of groups, using the same 
faint galaxy tracers as previous sections. 

To this aim, we define two samples of groups with at least 8 
galaxy members, one with groups located inside FVS and another 
with groups outside FVS. To make a direct comparison, we de- 
fine two new samples equating the mass distributions of the above 



Galaxy clustering 7 




in 
out 



I 



groups 
of g G M 




i i ■] in i | i i i 



10 12 14 

log(M) 



(a) 



10 12 14 

log(M) 



i.0.01 




Figure 4. (a) Mass distributions of galaxy groups. The left panel corresponds to the original samples of galaxy groups, where the dashed lines are the groups 
located inside FVS and the solid lines are the groups located outside FVS. The right panel corresponds to the new samples of groups redefined to have the 
same mass distributions. The bottom panels show the empirical cumulative mass distributions for each case, (b) Cross-correlation functions of galaxies in 
groups for gcM samples. The dashed lines correspond to galaxies in groups of gcM m sample and the solid lines correspond to galaxies in groups of gc,M ou , 
sample. Both samples have the same mass distributions. 




Figure 5. (a) The left panel shows the luminosity functions of the centre galaxies in groups, where these groups have the same mass distributions (i.e., 
luminosity distribution of samples gcM). The middle panel shows the luminosity distributions of the centre galaxies of the new samples redefined to have the 
same mass and luminosity distributions igoML). The right panel shows the mass distributions for groups composing samples ggML. This plot is to be sure 
that the mass distributions do not change when we adjusted the samples to have the same luminosity distributions. The bottom panels show the cumulative 
functions for each case. As the previous cases, in all panels the dashed lines indicate groups located in FVS, while the solid lines indicate groups located 
outside FVS. (b) Cross-correlation functions of galaxies in groups for gaML samples. The dashed lines correspond to galaxies in sample goML m and the 
solid lines correspond to galaxies in sample gaML,,,,,. Both samples have the same luminosity and mass distributions. 
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mentioned samples. The left panel of the Fig |6(a)| shows the mass 
distributions of groups with at least 8 members, where the dashed 
lines correspond to groups that are members of FVS and the solid 
lines correspond to groups that are not members of FVS. The right 
panel shows the mass distributions of the new samples, redefined 
to have the same mass distributions. The line coding is the same 
as in the previous sections. The bottom panels show the empirical 
cumulative mass distributions for each case. 

To determine the cross-correlation functions we use faint 
galaxies as tracers (defined in Section 14. Il l and the geometrical 
centres of groups as centres (G8M). Fig |6(b)| shows the cross- 
correlation functions for samples G8M. As in previous cases, the 
dashed lines correspond to centres of groups in sample G8M,„, 
while the solid lines correspond to centres of groups in sample 
GSM mll . As can be seen, there is no differences between the two 
correlation functions at scales of up to 2 tr 1 Mpc when virial 
masses are considered. 



4.7 Mock: g G ML - Mock and G8M - Mock samples. 

In order to assess the reproducibility of our previous results by cur- 
rent models for structure formation, we perform a similar analysis 
on a mock galaxy catalogue. We use the mock catalogue described 
in Section |2~4] which implements a semi-analytic model of galaxy 
formation in a ACDM cosmological framework. We have identi- 
fied both FVS and groups in this mock galaxy catalogue using the 
same procedures than that used in the SDSS-DR7. We estimate the 
cross-correlation functions for the corresponding samples to Sec- 
tion [43] and [476] The samples go ML - mock are selected from the 
mock catalogue following the same steps and taking into account 
the same restrictions implemented on the observational data (Sec- 
tion |4.5t . Fig. |7(a)| shows the corresponding cross-correlation func- 
tions, using the same line coding as in previous plots. Similarly, the 
samples G8M - mock are defined from the mock catalogue under 
the same conditions mentioned in Section |4~6l Fig. [7(b)] shows the 
cross-correlation functions obtained from these samples. As can be 
seen in these figures, the correlation functions estimated in mock 
catalogues are consistent with the results from observational data. 

4.8 Discussion 

The mass contained in the FVS is likely to form part of a virial- 
ized structure in the future, but can also be considered to be part 
of the original overdensity when it started its collapse in the past. 
Therefore, this mass can be used as a proxy for the equivalent peak 
height of the FVS. We will test this by comparing the amplitude of 
the correlation function at large separations, r > 10 Ir 1 Mpc, with 
that expected for peaks corresponding to the FVS masses. 

Fig. [8] shows the ratios between the correlation functions 
around objects inside and outside FVS, for G8M and G8M - mock 
samples. These ratios can be compared with the ratio between the 
corresponding bias factors expected for objects of the mass of the 
FVS and that of the centres outside superstructures. In order to do 
this, we need to obtain an estimate of the groups outside super- 
structures, and also of the FVSs. For the former, we use their virial 
masses, and for the latter, we measure their total luminosities and 
adopt the average mass-to-ligh t ratio for large-sc ale structures, 
< MIL >= 577 Mq/Lq, given bv lTinker et ail ( bold) . The bias fac- 
tors for the grou ps and FVSs are th en obtain ed using the fittin g 
formula given bv lTinker et all fcoioh . based on lsheth et all ( 1200 lh . 
The resulting ratio, ~ 2.4 ±0.3, is shown as a solid horizontal long- 



dashed line. As can be seen, in both real and mock data, these es- 
timates are consistent with the ratios between the measured cor- 
relation functions. Errors are represented by the horizontal dotted 
lines, located at the 16th and 84th percentils of the distribution of 
bias values for the set of haloes considered. 

The amplitude of the correlation function at large separations 
for objects lying inside FVS, is consistent with that expected for 
a structure of the mass of the FVS, and not with that expected for 
their virial mass. The interpretation for this result is that the po- 
sitions of the centres used for the measurement of the correlation 
function trace, on average, the location of the overdensity of the 
FVS. 

We test this hipothesis by measuring the centre of mass of 
the FVS using low and high mass groups living in them (we de- 
fine low mass groups as M < 10 13 M Q , and high mass groups as 
M > 10 14 M Q ). We find that the centres of mass are similar to 
within 1.1 IT 1 Mpc, which explains the small effect of the centres 
on the amplitude of the correlation function at large separations. 

We also measure the ratio between the correlation functions 
of these low and high mass groups living inside FVS. If they only 
traced the overdensity of the FVS, this ratio should be consistent 
with 1. Figure [9] shows this ratio in the top panel, confirming our 
hypothesis. The lower panel shows the same ratio taking as cen- 
tres all the low and high mass groups in the SDSS, and as can be 
seen, the ratio is significantly lower than 1, on average, as expected 
for this ratio in the case that the groups trace their own lower or 
higher overdensity corresponding to their lower or higher masses 
(respectively), rather than that of an FVS. 

As a final test, we also divided the sample of FVS in two, one 
for low and another for high FVS mass samples, and calculated the 
cross-correlation function between groups populating these sam- 
ples and tracer galaxies. As before, we ensure that the mass distri- 
butions of the groups are similar, to avoid the possible influence of 
the group mass on the amplitude of the correlation sample (even 
though we already showed that this is negligible when groups re- 
side in an FVS). We find indications that the clustering for groups 
in higher mass FVSs is higher, also lending more support to the 
results shown previously. 

Regarding smaller scales, those corresponding to the one-halo 
term, our analysis indicates that the differences in the clustering 
amplitude, taking into account whether galaxies reside in FVS or 
not, are the sum of several contributions. We identified at least two 
sources that may be responsible for these differences, i.e., the differ- 
ent clustering amplitude due to galaxy luminosity and galaxy host 
mass selection effects. We acknowledge that there may be other 
sources that contribute to this deviation related to the selection of 
distinct large-scale environments. However, our analysis indicates 
that deviations originated on those sources are not statistically sig- 
nificant, pointing at a scenario where the FVS environnment has no 
measurable effect on the 1-halo clustering. 



5 CONCLUSIONS 

In this work we present a detailed analysis of cross-correlation of 
faint and bright galaxies and its dependence on large-scale struc- 
tures. To this end, we use the properties of galax y groups to char- 
acterize local environments (Zapata et al. 2009), and th e proper- 
ties of future virialized structures dLuparello et alj|201lh as prox- 
ies of large scale environments. We distinguish between galaxies 
belonging to FVS and galaxies outside FVS in a volume limited 
sample up to z = 0.12. At large scales, the amplitude of the two 
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Figure 6. (a): The left panel shows the mass distributions of groups with at least 8 members. The right panel shows the mass distributions of the new samples 
of groups with at least 8 members, adjusted to have the same mass distributions (G8M). The bottom panels show the empirical cumulative mass distributions 
for each case. In all panels, the dashed lines correspond to groups inside FVS and the solid lines correspond to groups outside FVS. (b): Group-galaxy 
cross-correlation functions for samples G8M, where the groups have at least 8 members and have the same mass distributions. The dashed lines correspond 
to groups in sample G8M,„ and the solid lines correspond to groups in sample G8M„„ f . 
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Figure 7. (a): Mock: Cross-correlation functions of galaxies in groups. The dashed lines correspond to galaxies in sample gcML m - mock and the solid lines 
correspond to galaxies in sample gaML out — mock. Both samples have the same luminosity and mass distributions, (b): Mock: Group-galaxy cross-correlation 
functions, where the groups have at least 8 members and have the same mass distributions. The dashed lines correspond to groups in sample G8M,„ — mock 
and the solid lines correspond to groups in sample G8M OI „ — mock. 
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Figure 8. Ratio of the real-space correlation functions of galaxies inside 
and outside FVS, for samples G&M (upper panel) and GSM —mock (bottom 
panel), in solid lines, with uncertainties given by the shaded regions. The 
long-dashed horizontal lines shows the expected value for this ratio using 
the estimated FVS mass as that of the original overdensity that will eventu- 
ally give rise to a virialized structure, and that of the centres outside FVS. 
The dotted lines enclose the 1 - cr confidence values for this theoretical 
ratio. 
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Figure 9. Upper panel: cross-correlation function ratio, 
£in-low-mass(r)/f;in-hish-mass(r), for groups living inside FVS. Low 
mass corresponds to groups with M < 10 13 M Q while high mass corre- 
sponds to groups with M > 10 14 M G . Bottom panel: same as top panel, but 
for all the groups in SDSS separated in low and high mass samples. 



point cross-correlation function of bright and faint galaxies is sig- 
nificantly larger when the bright galaxies reside in FVS. At small 
scales, where the clustering cross-correlation signal is dominated 
by the local environment, the results inside and outside FVS are 
slightly different. In order to disentangle local and large scale con- 
tributions to the observed clustering, we have considered subsam- 
ples where the galaxy luminosity, and host group mass distributions 
are forced to be similar. Once these restrictions have been imple- 
mented, the resulting correlation functions show statistically negli- 
gible differences. 

In order to asses the reproducibility of these results within cur- 
rent models for structure formation, we have performed a similar 
analysis using a semi-analytic implementation in a ACDM cos- 
mological model. We have determined bright-faint galaxy cross- 
correlations taking into account the local and global environment 
using groups and FVS in a similar fashion as performed in the 
observations. By considering subsamples with similar galaxy lu- 
minosity and host group mass distributions we determine that the 
cross-correlations dependence on large-scale structures is consis- 
tent with the dependence found in the observations. 

This analysis suggests that the resulting behaviour of the 
cross-correlation clustering of bright and faint galaxies on large 
scale structures is a generic feature of galaxy clustering in hierar- 
chical scenarios and our current understanding of galaxy formation 
scenario. 

We studied the amplitude of the correlation function at large 
separation for centres inside and outside FVS, comparing them to 
that expected from the theory. We find that objects in FVS show 
a large-scale clustering consistent with that of an overdensity of 
the mass of the FVS in which they reside. These objects trace, on 
average, the centre of mass of the FVS, explaining this effect. Fur- 
thermore, for groups inside FVS, the amplitude of the correlation 
function at two-halo term separations does not seem to depend on 
group mass. This dependence is more clear when considering all 
the groups in the SDSS. 

In a forthcoming paper we will assess the properties of galax- 
ies in systems residing inside and outside FVS. 
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